A compendium of multi-omics data illuminating host responses to lethal human virus infections.

Human infections caused by viral pathogens trigger a complex gamut of host responses that limit disease, resolve infection, generate immunity, and contribute to severe disease or death. Here, we present experimental methods and multi-omics data capture approaches representing the global host response to infection generated from 45 individual experiments involving human viruses from the Orthomyxoviridae, Filoviridae, Flaviviridae, and Coronaviridae families. Analogous experimental designs were implemented across human or mouse host model systems, longitudinal samples were collected over defined time courses, and global multi-omics data (transcriptomics, proteomics, metabolomics, and lipidomics) were acquired by microarray, RNA sequencing, or mass spectrometry analyses. For comparison, we have included transcriptomics datasets from cells treated with type I and type II human interferon. Raw multi-omics data and metadata were deposited in public repositories, and we provide a central location linking the raw data with experimental metadata and ready-to-use, quality-controlled, statistically processed multi-omics datasets not previously available in any public repository. This compendium of infection-induced host response data for reuse will be useful for those endeavouring to understand viral disease pathophysiology and network biology.

Deciphering ApoE Genotype-Driven Proteomic and Lipidomic Alterations in Alzheimer's Disease Across Distinct Brain Regions.

Alzheimer's disease (AD) is a neurodegenerative disease with a complex etiology influenced by confounding factors such as genetic polymorphisms, age, sex, and race. Traditionally, AD research has not prioritized these influences, resulting in dramatically skewed cohorts such as three times the number of Apolipoprotein E (APOE) ε4-allele carriers in AD relative to healthy cohorts. Thus, the resulting molecular changes in AD have previously been complicated by the influence of apolipoprotein E disparities. To explore how apolipoprotein E polymorphism influences AD progression, 62 post-mortem patients consisting of 33 AD and 29 controls (Ctrl) were studied to balance the number of ε4-allele carriers and facilitate a molecular comparison of the apolipoprotein E genotype. Lipid and protein perturbations were assessed across AD diagnosed brains compared to Ctrl brains, ε4 allele carriers (APOE4+ for those carrying 1 or 2 ε4s and APOE4- for non-ε4 carriers), and differences in ε3ε3 and ε3ε4 Ctrl brains across two brain regions (frontal cortex (FCX) and cerebellum (CBM)). The region-specific influences of apolipoprotein E on AD mechanisms showcased mitochondrial dysfunction and cell proteostasis at the core of AD pathophysiology in the post-mortem brains, indicating these two processes may be influenced by genotypic differences and brain morphology.

PMart Web Application: Marketplace for Interactive Analysis of Panomics Data.

PMart is a web-based tool for reproducible quality control, exploratory data analysis, statistical analysis, and interactive visualization of 'omics data, based on the functionality of the pmartR R package. The newly improved user interface supports more 'omics data types, additional statistical capabilities, and enhanced options for creating downloadable graphics. PMart supports the analysis of label-free and isobaric-labeled (e.g., TMT, iTRAQ) proteomics, nuclear magnetic resonance (NMR) and mass-spectrometry (MS)-based metabolomics, MS-based lipidomics, and ribonucleic acid sequencing (RNA-seq) transcriptomics data. At the end of a PMart session, a report is available that summarizes the processing steps performed and includes the pmartR R package functions used to execute the data processing. In addition, built-in safeguards in the backend code prevent users from utilizing methods that are inappropriate based on omics data type. PMart is a user-friendly interface for conducting exploratory data analysis and statistical comparisons of omics data without programming.

Elucidating regulatory processes of intense physical activity by multi-omics analysis.

BACKGROUND: Physiological and biochemical processes across tissues of the body are regulated in response to the high demands of intense physical activity in several occupations, such as firefighting, law enforcement, military, and sports. A better understanding of such processes can ultimately help improve human performance and prevent illnesses in the work environment. METHODS: To study regulatory processes in intense physical activity simulating real-life conditions, we performed a multi-omics analysis of three biofluids (blood plasma, urine, and saliva) collected from 11 wildland firefighters before and after a 45 min, intense exercise regimen. Omics profiles post- versus pre-exercise were compared by Student's t-test followed by pathway analysis and comparison between the different omics modalities. RESULTS: Our multi-omics analysis identified and quantified 3835 proteins, 730 lipids and 182 metabolites combining the 3 different types of samples. The blood plasma analysis revealed signatures of tissue damage and acute repair response accompanied by enhanced carbon metabolism to meet energy demands. The urine analysis showed a strong, concomitant regulation of 6 out of 8 identified proteins from the renin-angiotensin system supporting increased excretion of catabolites, reabsorption of nutrients and maintenance of fluid balance. In saliva, we observed a decrease in 3 pro-inflammatory cytokines and an increase in 8 antimicrobial peptides. A systematic literature review identified 6 papers that support an altered susceptibility to respiratory infection. CONCLUSION: This study shows simultaneous regulatory signatures in biofluids indicative of homeostatic maintenance during intense physical activity with possible effects on increased infection susceptibility, suggesting that caution against respiratory diseases could benefit workers on highly physical demanding jobs.

malbacR: A Package for Standardized Implementation of Batch Correction Methods for Omics Data.

Mass spectrometry is a powerful tool for identifying and analyzing biomolecules such as metabolites and lipids in complex biological samples. Liquid chromatography and gas chromatography mass spectrometry studies quite commonly involve large numbers of samples, which can require significant time for sample preparation and analyses. To accommodate such studies, the samples are commonly split into batches. Inevitably, variations in sample handling, temperature fluctuation, imprecise timing, column degradation, and other factors result in systematic errors or biases of the measured abundances between the batches. Numerous methods are available via R packages to assist with batch correction for omics data; however, since these methods were developed by different research teams, the algorithms are available in separate R packages, each with different data input and output formats. We introduce the malbacR package, which consolidates 11 common batch effect correction methods for omics data into one place so users can easily implement and compare the following: pareto scaling, power scaling, range scaling, ComBat, EigenMS, NOMIS, RUV-random, QC-RLSC, WaveICA2.0, TIGER, and SERRF. The malbacR package standardizes data input and output formats across these batch correction methods. The package works in conjunction with the pmartR package, allowing users to seamlessly include the batch effect correction in a pmartR workflow without needing any additional data manipulation.

Multi-omics of NET formation and correlations with CNDP1, PSPB, and L-cystine levels in severe and mild COVID-19 infections.

The detailed mechanisms of COVID-19 infection pathology remain poorly understood. To improve our understanding of SARS-CoV-2 pathology, we performed a multi-omics and correlative analysis of an immunologically naïve SARS-CoV-2 clinical cohort from blood plasma of uninfected controls, mild, and severe infections. Consistent with previous observations, severe patient populations showed an elevation of pulmonary surfactant levels. Intriguingly, mild patients showed a statistically significant elevation in the carnosine dipeptidase modifying enzyme (CNDP1). Mild and severe patient populations showed a strong elevation in the metabolite L-cystine (oxidized form of the amino acid cysteine) and enzymes with roles in glutathione metabolism. Neutrophil extracellular traps (NETs) were observed in both mild and severe populations, and NET formation was higher in severe vs. mild samples. Our correlative analysis suggests a potential protective role for CNDP1 in suppressing PSPB release from the pulmonary space whereas NET formation correlates with increased PSPB levels and disease severity. In our discussion we put forward a possible model where NET formation drives pulmonary occlusions and CNDP1 promotes antioxidation, pleiotropic immune responses, and vasodilation by accelerating histamine synthesis.

Rapid remodeling of the soil lipidome in response to a drying-rewetting event.

BACKGROUND: Microbiomes contribute to multiple ecosystem services by transforming organic matter in the soil. Extreme shifts in the environment, such as drying-rewetting cycles during drought, can impact the microbial metabolism of organic matter by altering microbial physiology and function. These physiological responses are mediated in part by lipids that are responsible for regulating interactions between cells and the environment. Despite this critical role in regulating the microbial response to stress, little is known about microbial lipids and metabolites in the soil or how they influence phenotypes that are expressed under drying-rewetting cycles. To address this knowledge gap, we conducted a soil incubation experiment to simulate soil drying during a summer drought of an arid grassland, then measured the response of the soil lipidome and metabolome during the first 3 h after wet-up. RESULTS: Reduced nutrient access during soil drying incurred a replacement of membrane phospholipids, resulting in a diminished abundance of multiple phosphorus-rich membrane lipids. The hot and dry conditions increased the prevalence of sphingolipids and lipids containing long-chain polyunsaturated fatty acids, both of which are associated with heat and osmotic stress-mitigating properties in fungi. This novel finding suggests that lipids commonly present in eukaryotes such as fungi may play a significant role in supporting community resilience displayed by arid land soil microbiomes during drought. As early as 10 min after rewetting dry soil, distinct changes were observed in several lipids that had bacterial signatures including a rapid increase in the abundance of glycerophospholipids with saturated and short fatty acid chains, prototypical of bacterial membrane lipids. Polar metabolites including disaccharides, nucleic acids, organic acids, inositols, and amino acids also increased in abundance upon rewetting. This rapid metabolic reactivation and growth after rewetting coincided with an increase in the relative abundance of firmicutes, suggesting that members of this phylum were positively impacted by rewetting. CONCLUSIONS: Our study revealed specific changes in lipids and metabolites that are indicative of stress adaptation, substrate use, and cellular recovery during soil drying and subsequent rewetting. The drought-induced nutrient limitation was reflected in the lipidome and polar metabolome, both of which rapidly shifted (within hours) upon rewet. Reduced nutrient access in dry soil caused the replacement of glycerophospholipids with phosphorus-free lipids and impeded resource-expensive osmolyte accumulation. Elevated levels of ceramides and lipids with long-chain polyunsaturated fatty acids in dry soil suggest that lipids likely play an important role in the drought tolerance of microbial taxa capable of synthesizing these lipids. An increasing abundance of bacterial glycerophospholipids and triacylglycerols with fatty acids typical of bacteria and polar metabolites suggest a metabolic recovery in representative bacteria once the environmental conditions are conducive for growth. These results underscore the importance of the soil lipidome as a robust indicator of microbial community responses, especially at the short time scales of cell-environment reactions. Video Abstract.

pmartR 2.0: A Quality Control, Visualization, and Statistics Pipeline for Multiple Omics Datatypes.

The pmartR (https://github.com/pmartR/pmartR) package was designed for the quality control (QC) and analysis of mass spectrometry data, tailored to specific characteristics of proteomic (isobaric or labeled), metabolomic, and lipidomic data sets. Since its initial release, the tool has been expanded to address the needs of its growing userbase and now includes QC and statistics for nuclear magnetic resonance metabolomic data, and leverages the DESeq2, edgeR, and limma-voom R packages for transcriptomic data analyses. These improvements have made progress toward a unified omics processing pipeline for ease of reporting and streamlined statistical purposes. The package's statistics and visualization capabilities have also been expanded by adding support for paired data and by integrating pmartR with the trelliscopejs R package for the quick creation of trellis displays (https://github.com/hafen/trelliscopejs). Here, we present relevant examples of each of these enhancements to pmartR and highlight how each new feature benefits the omics community.

Simulated Night-Shift Schedule Disrupts the Plasma Lipidome and Reveals Early Markers of Cardiovascular Disease Risk.

Introduction: The circadian system coordinates daily rhythms in lipid metabolism, storage and utilization. Disruptions of internal circadian rhythms due to altered sleep/wake schedules, such as in night-shift work, have been implicated in increased risk of cardiovascular disease and metabolic disorders. To determine the impact of a night-shift schedule on the human blood plasma lipidome, an in-laboratory simulated shift work study was conducted. Methods: Fourteen healthy young adults were assigned to 3 days of either a simulated day or night-shift schedule, followed by a 24-h constant routine protocol with fixed environmental conditions, hourly isocaloric snacks, and constant wakefulness to investigate endogenous circadian rhythms. Blood plasma samples collected at 3-h intervals were subjected to untargeted lipidomics analysis. Results: More than 400 lipids were identified and quantified across 21 subclasses. Focusing on lipids with low between-subject variation per shift condition, alterations in the circulating plasma lipidome revealed generally increased mean triglyceride levels and decreased mean phospholipid levels after night-shift relative to day-shift. The circadian rhythms of triglycerides containing odd chain fatty acids peaked earlier during constant routine after night-shift. Regardless of shift condition, triglycerides tended to either peak or be depleted at 16:30 h, with chain-specific differences associated with the direction of change. Discussion: The simulated night-shift schedule was associated with altered temporal patterns in the lipidome. This may be premorbid to the elevated cardiovascular risk that has been found epidemiologically in night-shift workers.

Hypergraph models of biological networks to identify genes critical to pathogenic viral response.

BACKGROUND: Representing biological networks as graphs is a powerful approach to reveal underlying patterns, signatures, and critical components from high-throughput biomolecular data. However, graphs do not natively capture the multi-way relationships present among genes and proteins in biological systems. Hypergraphs are generalizations of graphs that naturally model multi-way relationships and have shown promise in modeling systems such as protein complexes and metabolic reactions. In this paper we seek to understand how hypergraphs can more faithfully identify, and potentially predict, important genes based on complex relationships inferred from genomic expression data sets. RESULTS: We compiled a novel data set of transcriptional host response to pathogenic viral infections and formulated relationships between genes as a hypergraph where hyperedges represent significantly perturbed genes, and vertices represent individual biological samples with specific experimental conditions. We find that hypergraph betweenness centrality is a superior method for identification of genes important to viral response when compared with graph centrality. CONCLUSIONS: Our results demonstrate the utility of using hypergraphs to represent complex biological systems and highlight central important responses in common to a variety of highly pathogenic viruses.

A resource of lipidomics and metabolomics data from individuals with undiagnosed diseases.

Every year individuals experience symptoms that remain undiagnosed by healthcare providers. In the United States, these rare diseases are defined as a condition that affects fewer than 200,000 individuals. However, there are an estimated 7000 rare diseases, and there are an estimated 25-30 million Americans in total (7.6-9.2% of the population as of 2018) affected by such disorders. The NIH Common Fund Undiagnosed Diseases Network (UDN) seeks to provide diagnoses for individuals with undiagnosed disease. Mass spectrometry-based metabolomics and lipidomics analyses could advance the collective understanding of individual symptoms and advance diagnoses for individuals with heretofore undiagnosed disease. Here, we report the mass spectrometry-based metabolomics and lipidomics analyses of blood plasma, urine, and cerebrospinal fluid from 148 patients within the UDN and their families, as well as from a reference population of over 100 individuals with no known metabolic diseases. The raw and processed data are available to the research community so that they might be useful in the diagnoses of current or future patients suffering from undiagnosed disorders.

Proteogenomic and metabolomic characterization of human glioblastoma.

Glioblastoma (GBM) is the most aggressive nervous system cancer. Understanding its molecular pathogenesis is crucial to improving diagnosis and treatment. Integrated analysis of genomic, proteomic, post-translational modification and metabolomic data on 99 treatment-naive GBMs provides insights to GBM biology. We identify key phosphorylation events (e.g., phosphorylated PTPN11 and PLCG1) as potential switches mediating oncogenic pathway activation, as well as potential targets for EGFR-, TP53-, and RB1-altered tumors. Immune subtypes with distinct immune cell types are discovered using bulk omics methodologies, validated by snRNA-seq, and correlated with specific expression and histone acetylation patterns. Histone H2B acetylation in classical-like and immune-low GBM is driven largely by BRDs, CREBBP, and EP300. Integrated metabolomic and proteomic data identify specific lipid distributions across subtypes and distinct global metabolic changes in IDH-mutated tumors. This work highlights biological relationships that could contribute to stratification of GBM patients for more effective treatment.

Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR-MS) peak intensity normalization for complex mixture analyses.

RATIONALE: Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR-MS) is a preferred technique for analyzing complex organic mixtures. Currently, there is no consensus normalization approach, nor an objective method for selecting one, for quantitative analyses of FT-ICR-MS data. We investigate a method to evaluate and score the amount of bias various normalization approaches introduce into the data. METHODS: We evaluate the ability of the Statistical Procedure for the Analysis of Normalization Strategies (SPANS) to guide the selection of appropriate normalization approaches for two different FT-ICR-MS data sets. Furthermore, we test the robustness of SPANS results to changes in SPANS parameter values and assess the impact of using various normalization approaches on downstream statistical analyses. RESULTS: The normalization approach identified by SPANS differed for the two data sets. Normalization methods impacted the statistical significance of peaks differently, underscoring the importance of carefully evaluating potential methods. More consistent SPANS scores resulted when at least 120 significant peaks are used, where larger sets of peaks were obtained by increasing the p-value threshold. Interestingly, we show that total sum scaling and highest peak normalization, used in previous studies, underperformed relative to SPANS-recommended normalization approaches. CONCLUSIONS: Although there is no single, best normalization method for all data sets, SPANS provides a mechanism to identify an appropriate normalization method for analyzing FT-ICR-MS data quantitatively. The number of peaks used in the background distributions of SPANS contributes more significantly to the reproducibility of results than the p-value thresholds used to obtain those peaks.

Correction: yippee like 3 (ypel3) is a novel gene required for myelinating and perineurial glia development.

[This corrects the article DOI: 10.1371/journal.pgen.1008841.].

Unveiling molecular signatures of preeclampsia and gestational diabetes mellitus with multi-omics and innovative cheminformatics visualization tools.

To fully enable the development of diagnostic tools and progressive pharmaceutical drugs, it is imperative to understand the molecular changes occurring before and during disease onset and progression. Systems biology assessments utilizing multi-omic analyses (e.g. the combination of proteomics, lipidomics, genomics, etc.) have shown enormous value in determining molecules prevalent in diseases and their associated mechanisms. Herein, we utilized multi-omic evaluations, multi-dimensional analysis methods, and new cheminformatics-based visualization tools to provide an in depth understanding of the molecular changes taking place in preeclampsia (PRE) and gestational diabetes mellitus (GDM) patients. Since PRE and GDM are two prevalent pregnancy complications that result in adverse health effects for both the mother and fetus during pregnancy and later in life, a better understanding of each is essential. The multi-omic evaluations performed here provide new insight into the end-stage molecular profiles of each disease, thereby supplying information potentially crucial for earlier diagnosis and treatments.

A global lipid map defines a network essential for Zika virus replication.

Zika virus (ZIKV), an arbovirus of global concern, remodels intracellular membranes to form replication sites. How ZIKV dysregulates lipid networks to allow this, and consequences for disease, is poorly understood. Here, we perform comprehensive lipidomics to create a lipid network map during ZIKV infection. We find that ZIKV significantly alters host lipid composition, with the most striking changes seen within subclasses of sphingolipids. Ectopic expression of ZIKV NS4B protein results in similar changes, demonstrating a role for NS4B in modulating sphingolipid pathways. Disruption of sphingolipid biosynthesis in various cell types, including human neural progenitor cells, blocks ZIKV infection. Additionally, the sphingolipid ceramide redistributes to ZIKV replication sites, and increasing ceramide levels by multiple pathways sensitizes cells to ZIKV infection. Thus, we identify a sphingolipid metabolic network with a critical role in ZIKV replication and show that ceramide flux is a key mediator of ZIKV infection.

yippee like 3 (ypel3) is a novel gene required for myelinating and perineurial glia development.

Hypomyelination, a neurological condition characterized by decreased production of myelin sheets by glial cells, often has no known etiology. Elucidating the genetic causes of hypomyelination provides a better understanding of myelination, as well as means to diagnose, council, and treat patients. Here, we present evidence that YIPPEE LIKE 3 (YPEL3), a gene whose developmental role was previously unknown, is required for central and peripheral glial cell development. We identified a child with a constellation of clinical features including cerebral hypomyelination, abnormal peripheral nerve conduction, hypotonia, areflexia, and hypertrophic peripheral nerves. Exome and genome sequencing revealed a de novo mutation that creates a frameshift in the open reading frame of YPEL3, leading to an early stop codon. We used zebrafish as a model system to validate that YPEL3 mutations are causative of neuropathy. We found that ypel3 is expressed in the zebrafish central and peripheral nervous system. Using CRISPR/Cas9 technology, we created zebrafish mutants carrying a genomic lesion similar to that of the patient. Our analysis revealed that Ypel3 is required for development of oligodendrocyte precursor cells, timely exit of the perineurial glial precursors from the central nervous system (CNS), formation of the perineurium, and Schwann cell maturation. Consistent with these observations, zebrafish ypel3 mutants have metabolomic signatures characteristic of oligodendrocyte and Schwann cell differentiation defects, show decreased levels of Myelin basic protein in the central and peripheral nervous system, and develop defasciculated peripheral nerves. Locomotion defects were observed in adult zebrafish ypel3 mutants. These studies demonstrate that Ypel3 is a novel gene required for perineurial cell development and glial myelination.

ftmsRanalysis: An R package for exploratory data analysis and interactive visualization of FT-MS data.

The high-resolution and mass accuracy of Fourier transform mass spectrometry (FT-MS) has made it an increasingly popular technique for discerning the composition of soil, plant and aquatic samples containing complex mixtures of proteins, carbohydrates, lipids, lignins, hydrocarbons, phytochemicals and other compounds. Thus, there is a growing demand for informatics tools to analyze FT-MS data that will aid investigators seeking to understand the availability of carbon compounds to biotic and abiotic oxidation and to compare fundamental chemical properties of complex samples across groups. We present ftmsRanalysis, an R package which provides an extensive collection of data formatting and processing, filtering, visualization, and sample and group comparison functionalities. The package provides a suite of plotting methods and enables expedient, flexible and interactive visualization of complex datasets through functions which link to a powerful and interactive visualization user interface, Trelliscope. Example analysis using FT-MS data from a soil microbiology study demonstrates the core functionality of the package and highlights the capabilities for producing interactive visualizations.

Automated mass spectrometry imaging of over 2000 proteins from tissue sections at 100-µm spatial resolution.

Biological tissues exhibit complex spatial heterogeneity that directs the functions of multicellular organisms. Quantifying protein expression is essential for elucidating processes within complex biological assemblies. Imaging mass spectrometry (IMS) is a powerful emerging tool for mapping the spatial distribution of metabolites and lipids across tissue surfaces, but technical challenges have limited the application of IMS to the analysis of proteomes. Methods for probing the spatial distribution of the proteome have generally relied on the use of labels and/or antibodies, which limits multiplexing and requires a priori knowledge of protein targets. Past efforts to make spatially resolved proteome measurements across tissues have had limited spatial resolution and proteome coverage and have relied on manual workflows. Here, we demonstrate an automated approach to imaging that utilizes label-free nanoproteomics to analyze tissue voxels, generating quantitative cell-type-specific images for >2000 proteins with 100-µm spatial resolution across mouse uterine tissue sections preparing for blastocyst implantation.

Corrigendum: Light-Stress Influences the Composition of the Murine Gut Microbiome, Memory Function, and Plasma Metabolome.

[This corrects the article DOI: 10.3389/fmolb.2019.00108.].